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While techniques to compute thermal fluctuation induced, or pseudo-Casimir, forces in equilibrium 
systems are well established, the same is not true for non-equilibrium cases. We present a general 
formalism that allows us to unambiguously compute non-equilibrium fluctuation induced forces by 
specifying the energy of interaction of the fluctuating fields with the boundaries. For a general class 
of classical fields with dissipative dynamics, we derive a very general relation between the Laplace 
transform of the time-dependent force and the static partition function for a related problem with 
a different Hamiltonian. In particular, we demonstrate the power of our approach by computing, 
for the first time, the explicit time dependence of the non-equilibrium pseudo-Casimir force induced 
between two parallel plates, upon a sudden change in the temperature of the system. We also show 
how our results can be used to determine the steady-state behavior of the non-equilibrium force in 
systems where the fluctuations are driven by colored noise. 
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The Casimir effect arises when the fluctuations of a 
quantum or classical field are constrained by the presence 
of surfaces or objects placed in the field IH |2l |3] . While 
the standard Casimir effect refers to interactions arising 
from constrained fluctuations of an electromagnetic field, 
a similar effect that arises due to constrained thermal 
fluctuations in systems with long range correlations, such 
as critical fluids, smectic manifolds or liquid crystals, is 
termed the pseudo-Casimir effect. The simplest system 
where the so called pseudo-Casimir effect arises is the 
classical free scalar field theory, and long range pseudo- 
Casimir forces arise when the field theory is massless, i.e. 
where the Hamiltonian of the system is given by: 

H=\jd^ [Vmf- (1) 

The fluctuating fleld (f> can describe the order parameter 
for critical systems, such as binary liquids at the critical 
point [4J , and the occurrence of a pseudo-Casimir interac- 
tion in such a system has recently been confirmed exper- 
imentally [5] . The field (/) can also represent the phase of 
the complex order parameter in a superfluid state, where 
pseudo-Casimir forces can lead to the thinning of super- 
fluid ^He fllms [6]. Free vectorial field theories describing 
liquid crystal systems also exhibit pseudo-Casimir type 
interactions between the surfaces confining the system 
[7]. The interaction induced between surfaces or objects 
immersed in such systems can be considered to be due to 
the imposition of boundary conditions on the field or due 
to an energy of interaction of the surfaces or objects with 
the field. If, for instance, we consider the surfaces to be 
two parallel plates and specify the boundary conditions 
for the field on the two plates {e.g. Dirichlet, Neumann 
or Robin), the pseudo-Casimir interaction between the 
plates can be computed from the free energy using stan- 
dard techniques if the system is at thermal equilibrium 



[T]. The same techniques cannot, however, be used to 
compute the pseudo-Casimir force if the system is driven 
out of equilibrium, for example by a sudden change in 
temperature or by colored noise forcing. Since many ex- 
perimental and naturally occurring systems, where such 
fiuctuation induced forces are important {e.g. between 
inclusions in cell membranes) , are actually out of equilib- 
rium, it is essential to have a well-defined, unambiguous 
way to compute non-equilibrium pseudo-Casimir forces. 
One of the principal problems when analyzing the non- 
equilibrium pseudo-Casimir effect, say, for parallel plates, 
is to obtain a correct expression for the force between the 
two plates. In previous studies, the stress tensor has been 
used to study both the dynamical behavior of the force 
O [9l [To] and the force fluctuations in equilibrium 
However, only the average force in equilibrium can be 
strictly computed using the stress tensor. Results using 
the stress tensor may, however, be reliable for situations 
close to equilibrium [TOj . Another approach is to con- 
struct a model with a specified non equilibrium dynamics 
and to specify by hand the force at the wall. For exam- 
ple, in [T5], the dynamical field was related to a particle 
density and the local pressure on the wall is then given by 
the ideal gas form via kinetic reasoning. In this letter, 
rather than imposing boundary conditions on the field 
we specify the energy of interaction of the wall with the 
field. In this way we can write down the instantaneous 
force on the wall unambiguously. Our novel approach 
can be used to recover results for the usual boundary 
conditions employed in studies of the pseudo-Casimir in- 
teraction by taking the appropriate limit but is, in fact, 
extremely general. We believe that such a microscopic 
approach is required to obtain meaningful results for the 
pseudo-Casimir force out of equilibrium. 

We commence by considering the most general case of 



2 



a free field theory where the Hamiltonian can be written 
in terms of a general quadratic Hamiltonian. 



H = 



1 



dxdx'(/)(x)A(x,x',/)0(x'), 



(2) 



where A is aself-adoint operator i.e. A(x,x') = A(x',x). 
Here I represents any suitable free parameter in the prob- 
lem but, for concreteness, it could be the separation of 
two parallel plates which interact with the field. Now, if 
we choose 

A(x, x', I) = [-V^ + S{z)c + S{z - l)c\ (5(x - x'), (3) 

then this corresponds to a free field theory where the 
fluctuations of the field (j) are suppressed at both plates 
(at z — and ai z = I). Clearly, when c — s- oo one will 
obtain Dirichlet boundary conditions at the two plates. 
The instantaneous generalized force acting on the plate 
at z = Z is then given by 



dH 

'In 



^ I dxdx'0(x)|A(x,x',O0(x'). 



(4) 



The equilibrium value of this force can also be written in 
the familiar form 



{F,)=T^^HZ{A)) 



(5) 



where Z{A) = J d[(f)] exp(— [A]) with H as defined in 
Eq.(|2]) and where /3 = l/T is the inverse temperature. 
Alternatively the force can be expressed via Eq.(|4| as 



(Fi) 



T 

— / dxdx' 



|a(x,x',0 



A-i(x,x',0 (6) 



We now consider the dynamical problem where the sys- 
tem is prepared in a state (j) — at the time t — (this 
could have been by cooling the system to a very low tem- 
perature for instance) and then letting it relax at some 
non-zero temperature T. We will consider the very gen- 
eral relaxational dynamics, where the evolution of the 
field is given by 



dt 



-i?A0(x,t) -K7?(x,t) 



(7) 



where the noise is Gaussian and uncorrelated in time 
mean with correlation function {ri(Ti^t)ri(x' ,t')) — 
2T6{t — t')R{x — x'). Here i? is a symmetric transla- 
tionally invariant operator and the choice of the noise 
correlator ensures that the dynamics obeys detailed bal- 
ance. For thermal noise uncorrelated in time this is a 
very general representation of the dynamics of soft con- 
densed matter systems, such as those mentioned above, 
when intertial and relativistic effects (due to the pres- 
ence of large viscosities) can be neglected. The formal 
solution to this equation, for our initial condition, is 



(8) 



x,t) = / dsexp (— (i — s)i?A) 77(x, s) 
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FIG. 1: Approach to equilibrium for various values of the 
plate separation I (top to bottom I — 1,2,4) for d — 3. 
Nonequilibrium force F, obtained by the direct numerical in- 



verse Laplace transform of Eq. (14 1, plotted in units of the 
equilibrium pseudo-Casimir force (F^q) as a function of time. 

This means that the equal time correlation function of 
the field, C(x,x',i) = {4>{x,t)(f>{x' ,t)), can be shown to 
be given by 



C{t) =TR[1- exp(-2iAi?)] (AR)- 



(9) 



Now if we Laplace transform this equation (defining 
£/(s) — dtexp{—st)f{t))we find that 



/:C(x,x',s) = -A7i(x,x,n 
s 



where 



A,, = A - 



sR- 



(10) 



(11) 



Now using this and Eq.([4| we find that the Laplace 
transform of the average value of the generalised force 
is given by 



{CFi{s))^~^^ J dxdx'|A(x,x',/)A;i(x,x'), (12) 
which using Eq, ([TTJ and Eq. ^) can be written as 
(CFiis)) = -§/ dxdx'|A,(x,x',OA;i(x,x') 



T d 
75/ 



ln(Z(A,)) 



(13) 



This result is quite remarkable. It means that the Laplace 
transform of the time dependent pseudo-Casimir force 
considered here is given by a static pseudo-Casimir force 
with the operator A^. It is clear that the result is also 
valid for the force on any surface in the system which 
interacts with the field. Providing the static partition 
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function is known for the screened problem, the corre- 
sponding time dependent force can be extracted by in- 
verting the Laplace transform. 

We now turn to the case where the imposed bound- 
ary conditions are Dirichlet and the two plates are im- 
mersed in the fluctuating medium (hence there is fluctu- 
ating medium on both sides of each plate) . For the relax- 
ational dynamics of the field (Eqj?]) we will consider the 
simplest case with R= (5(x — x'), which means that the 
Gaussian noise is both uncorrelated in space and time. 
This corresponds to dynamics that does not conserve the 
order parameter and is also known as model A dynamics. 
We take the total length of the system to be L which is 
fixed, and place the plates of area A a distance I apart. 
Standard results on the screened pseudo-Casimir inter- 
action [T] give 

2AT 

where d is the dimension of the space, A is the area of 
the plates and. 



"exp(-2?yfc2T|) 



1 



(15) 



The equilibrium behavior is easily extracted by examin- 
ing the pole at s = which yields, as anticipated, the 
standard equilibrium pseudo-Casimir force [13] 

ATT{d)C{d) 



(16) 



where F is Euler's gamma function and C is the Rie- 
mann zeta function [14 . The full time dependence of the 
force, starting at zero at t = and relaxing to the equi- 
librium value above, can be extracted by direct Laplace 
inversion of Eq. (14). Fig.Q shows the approach to 
equilibrium for three different plate separations. Clearly, 
the relaxation times increase with plate separation and 
this is due to the fact that the underlying dynamics is 
diffusive and hence P sets a time scale. Useful analytic 
expressions for the early and late time behavior of the 



non-equilibrium force can also be obtained from Eq.(14l 



Using specific properties of Laplace transforms and the 
fact that Fi{t = 0) = allows us to arrive at an expres- 
sion for the time derivative of the force 



' dt ' 



AT 



,0 -.d d + 2 



E 



t 



1 



exp(- 



which in turn gives us the asymptotic behavior of the 
pseudo-Casimir force for different time regimes defined 
by P/t 

im)) - -7^^exp(--) for - »1 



AT 



d(87r)*t5 



P 

for — < 1 



(18) 



It is interesting to note that the late time correction is in- 
dependent of I. This is because the medium between the 
two plates has a relaxation time t{1) ^ P j^iP whereas 
the slowest relaxation times in the system are associated 
with the medium outside the two plates and hence at late 
times the correction is dominated by the relaxation of 
the external system in the thermodynamic limit L ^ oo. 
This diffusive relaxation is responsible for the power law 
approach to equilibrium. It is natural to ask how these 
model A results translate into the stress tensor formal- 
ism. Because the rhs of Eq. ( 13 1 is static we can write 



the result in terms of the stress tensor for the correspond- 
ing field theory, and inverting the Laplace transform we 
find |,16] an effective dynamical stress tensor: 



T 



dyn 



1 



l_d_ 
2di' 



(19) 



Thus for the problem considered here we see that the 
first part of the stress tensor picks up a time derivative 
term whose expectation will vanish in equilibrium to give 
the usual equilibrium stress tensor result. It is also clear 
that simply using the standard stress tensor on a surface 
in a system out of equilibrium will only give the right 
result when Dirichlet boundary conditions arc imposed 
(as = on the surface). 

The method can also be used to examine the dynamics 
resulting from a sudden change in temperature, from say 
To where the system is in equilibrium to a temperature T. 
In this case, the initial configuration of the field (/)(x, 0) 
has the correlation function 



(<^(x,O)0(x',O)) = roA"'(x,x',O. 



(20) 



Solving the equation of motion Eq.Q with this initial 
condition yields the time dependent correlation function 

C(x, x', t) ^ To cxp(-2tA)A"i + T{1 - exp(-2tA))A"i 

(21) 

The same reasoning as for Eq.([9|-(13l yields 



T-n d 
~s di 



ln(Z(A,)) (22) 



where we have used the fact that ^ln(Z(As)) is inde- 
pendent of the temperature. For Dirichlet boundary con- 



ditions this gives, in analogy with Eq.(18l, the limiting 
behavior 



{Flit)) ~ {Fi),gTo- 



2A{T-To) , P. , P , 
— ^ ^exp(- — ) for — >1 



P 



d[SiT)^t2 t 



(23) 



One can also consider the behavior of the pseudo- 
Casimir force for relaxational dynamics of the form of 
Eq.(|7]) but where the forcing noise is colored in time such 
that (r/(x,t)r/(x',i')) = T5(x - x')u; exp(-u;|t - t'|) (so 
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FIG. 2: Nonequilibrium pseudo-Casimir force, F, in units of 
the equilibrium force (Feg) as a function of uj for / — 1 and 
rf = 3. Inset: F, as a function of plate separation / for lo = 1 
and d = 3 (solid line). The equilibrium force (dashed line) is 
shown for comparison. 

i?(x — x') = (5(x — x') has the model A from). Here T 
represents an energy scale, u a frequency and the result- 
ing steady state is not an equilibrium one. The average 
value of the force in the steady state regime can be com- 
puted using the same formalism above and we find that 
the correlation function of the field is given by 

C(x,x',c^) = T [A(x,x',0-^ - A2^(x,x',0-^] (24) 

which yields 



(Fiico)) = T- [ln(Z(A)) - HZ{A,^)] (25) 

Hence again we find that one can compute a force in a 
non equilibrium system from knowledge of static screened 
systems [T^. Note that in the limit w ^ oo we recover 
the white noise equilibrium result of Eq. ([s]). Fig. ^ 
shows the frequency dependence of the non-equilibrium 



force obtained from Eq. ( 25 1 for a two plate system with 
Dirichlet boundary conditions as we had before. Again 
P sets a timescale and we see that for uj the force, 

F tends to the equilibrium white noise value, F^q, as ex- 
pected, while for uj <C l^^, F <C F^q and as a; — ^ 0, F 
vanishes. The inset to Fig. ^ shows how the force de- 
pends on plate separation for fixed uj. Again, equilibrium 
behavior is recovered for large I {uj ^ while for small 
plate separations the force changes qualitatively scaling 
as We note that the result Eq.(25l agrees with a 



computation for the same system where the steady state 
force was computed using the stress tensor . We finally 
note that our approach can also be readily extended to 
study the non-equilibrium pseudo-Casimir force between 
two small defect regions within the pairwise approxima- 
tion [16,. 

Previous studies on the dynamical pseudo-Casimir ef- 
fect concentrated on steady state non equilibrium dy- 
namics or dynamics close to equilibrium and considered 



Dirichlet boundary conditions assuming that the equi- 
librium stress tensor could be applied to compute the 
force. Our formalism marks a major advance that over- 
comes these restrictions by allowing the time-dependent 
force to be evaluated unambiguously via an expression 
for the energy of the field. We emphasize the general- 
ity of our approach as it is valid for (i) any dissipative 
dynamics of the form Eq. ([t]) , (ii) any free field the- 
ory {e.g. operators A which contain terms such as 
as is the case for the height fluctuations of lipid mem- 
branes) and (iii) for any generalized force conjugate to 
any parameter in the system. It is worth noticing that 
the static theory that must be solved for conserved dy- 
namics (model B), where R = — V^, becomes non-local 
and the corresponding static calculation presents us with 
the interesting problem of understanding pseudo-Casimir 
forces for systems with non local interactions. This re- 
search was supported in part by the National Science 
Foundation under Grant No.PHY05-5II64 (while at the 
KITP UCSB program The theory and practice of fluc- 
tuation induced interactions 2008). DSD acknowledges 
support from the Institut Universitaire de France. 



[1] 
[2] 
[3] 

[4 
[5 

[6: 

[7 

[s: 

[9 
[10 

[11 
[12 

[13 
[14 

[15 
[16 



M. Kardar and R. Golestanian, Rev. Mod. Phys. 71, 1233 
(1999) 

V.M. Mostepanenko and N.N. Trunov, The Casimir Ef- 
fect and its Applications, (Oxford) (1997) 
K. A. Milton, The Casimir Effect: Physical Manifesta- 
tions of Zero-Point Energy (World Scientific, Singapore) 
(2001). 

M.E. Fisher and P.-G. de Gennes, C. R. Acad. Sci. Paris 
B 287, 207 (1978) 

C. Hertlein, L. Helden, A. Gambassi, S. Dietrich and C. 
Bechinger Nature 451, 172 (2008) 

R. Zandi, J. Rudnick and M. Kardar, Phys. Rev. Lett. 
93, 155302 (2004) 

A. Ajdari, L. Peliti, and J. Prost, Phys. Rev. Lett. 66, 
1481 (1991) 

D. Bartolo, A. Adjari and J.-B. Fournier, Phys. Rev. E 
67, 061112 (2003) 

A. Najafi and R. Golestanian, Europhys. Lett., 68 , 776 
(2004) 

A. Gambassi and S. Dietrich, J. Stat. Phys. 123, 929 
(2006) 

D. Bartolo, A. Adjari, J.-B. Fournier and R. Golestanian, 
Phys. Rev. Lett. 89, 230601 (2002) 

R. Brito, U. Marini Bettolo Marconi and R. Soto, Phys. 
Rev E, 76, 011113 (2007) 

M. Krech and S. Dietrich, Phys. Rev. A 46, 1886 (1992) 
I.S. Gradshteyn and I.M. Ryzhik (A. Jeffrey and D. Zwill- 
inger (eds)). Table of Integrals Series and Products, (Aca- 
demic Press) (2000) 

This link between the screened equilibrium Casimir force 
and that arising due to colored noise forcing for parallel 
plates with Dirichlet boundary was first noticed in [8], 
we show here that it is a very general phenomenon. 
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